Human Endogenous Retrovirus (HERV) Transcriptome Is Dynamically Modulated during SARS-CoV-2 Infection and Allows Discrimination of COVID-19 Clinical Stages.

SARS-CoV-2 infection is known to trigger an important inflammatory response, which has a major role in COVID-19 pathogenesis. In infectious and inflammatory contexts, the modulation of human endogenous retroviruses (HERV) has been broadly reported, being able to further sustain innate immune responses due to the expression of immunogenic viral transcripts, including double-stranded DNA (dsRNA), and eventually, immunogenic proteins. To gain insights on this poorly characterized interplay, we performed a high-throughput expression analysis of ~3,300 specific HERV loci in the peripheral blood mononuclear cells (PBMCs) of 10 healthy controls and 16 individuals being either convalescent after the infection (6) or retesting positive after convalescence (10) (Gene Expression Onmibus [GEO] data set GSE166253). Results showed that the exposure to SARS-CoV-2 infection modulates HERV expression according to the disease stage and reflecting COVID-19 immune signatures. The differential expression analysis between healthy control (HC) and COVID-19 patients allowed us to identify a total of 282 differentially expressed HERV loci (deHERV) in the individuals exposed to SARS-CoV-2 infection, independently from the clinical form. In addition, 278 and 60 deHERV loci that were specifically modulated in individuals convalescent after COVID19 infection (C) and patients that retested positive to SARS-CoV-2 after convalescence (RTP) as individually compared to HC, respectively, as well as 164 deHERV loci between C and RTP patients were identified. The identified HERV loci belonged to 36 different HERV groups, including members of all three classes. The present study provides an exhaustive picture of the HERV transcriptome in PBMCs and its dynamic variation in the presence of COVID-19, revealing specific modulation patterns according to the infection stage that can be relevant to the disease clinical manifestation and outcome. IMPORTANCE We report here the first high-throughput analysis of HERV loci expression along SARS-CoV-2 infection, as performed with peripheral blood mononuclear cells (PBMCs). Such cells are not directly infected by the virus but have a crucial role in the plethora of inflammatory and immune events that constitute a major hallmark of COVID-19 pathogenesis. Results provide a novel and exhaustive picture of HERV expression in PBMCs, revealing specific modulation patterns according to the disease condition and the concomitant immune activation. To our knowledge, this is the first set of deHERVs whose expression is dynamically modulated across COVID-19 stages, confirming a tight interplay between HERV and cellular immunity and revealing specific transcriptional signatures that can have an impact on the disease clinical manifestation and outcome.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans.

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV's roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species' genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates' germlines after the segregation of Gibbons' species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold-fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells' populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein-protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer's, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells' differentiation programs.

Human endogenous retrovirus K in the respiratory tract is associated with COVID-19 physiopathology.

BACKGROUND: Critically ill 2019 coronavirus disease (COVID-19) patients under invasive mechanical ventilation (IMV) are 10 to 40 times more likely to die than the general population. Although progression from mild to severe COVID-19 has been associated with hypoxia, uncontrolled inflammation, and coagulopathy, the mechanisms involved in the progression to severity are poorly understood. METHODS: The virome of tracheal aspirates (TA) from 25 COVID-19 patients under IMV was assessed through unbiased RNA sequencing (RNA-seq), and correlation analyses were conducted using available clinical data. Unbiased sequences from nasopharyngeal swabs (NS) from mild cases and TA from non-COVID patients were included in our study for further comparisons. RESULTS: We found higher levels and differential expression of human endogenous retrovirus K (HERV-K) genes in TA from critically ill and deceased patients when comparing nasopharyngeal swabs from mild cases to TA from non-COVID patients. In critically ill patients, higher HERV-K levels were associated with early mortality (within 14 days of diagnosis) in the intensive care unit. Increased HERV-K expression in deceased patients was associated with IL-17-related inflammation, monocyte activation, and an increased consumption of clotting/fibrinolysis factors. Moreover, increased HERV-K expression was detected in human primary monocytes from healthy donors after experimental SARS-CoV-2 infection in vitro. CONCLUSION: Our data implicate the levels of HERV-K transcripts in the physiopathology of COVID-19 in the respiratory tract of patients under invasive mechanical ventilation. Video abstract.

High Expression of HERV-K (HML-2) Might Stimulate Interferon in COVID-19 Patients.

Background. Interferon is a marker of host antiviral immunity, which is disordered in COVID-19 patients. ERV can affect the secretion of interferon through the cGAS-STING pathway. In this study, we explored whether IFN-I and HERV-K (HML-2) were activated in COVID-19 patients and whether there was an interaction between them. Methods. We collected blood samples from COVID-19 patients and healthy controls. We first detected the expression of HERV-K (HML-2) gag, env, and pol genes and IFN-I-related genes between patients and healthy people by qPCR, synchronously detected VERO cells infected with SARS-CoV-2. Then, the chromosome distributions of highly expressed HERV-K (HML-2) gag, env, and pol genes were mapped by the next-generation sequencing results, and GO analysis was performed on the related genes. Results. We found that the HERV-K (HML-2) gag, env, and pol genes were highly expressed in COVID-19 patients and VERO cells infected with SARS-CoV-2. The interferon-related genes IFNB1, ISG15, and IFIT1 were also activated in COVID-19 patients, and GO analysis showed that HERV-K (HML-2) can regulate the secretion of interferon. Conclusions. The high expression of HERV-K (HML-2) might activate the increase of interferon in COVID-19 patients, proving that HERV-K does not only play a negative role in the human body.

Rare isolation of human-tropic recombinant porcine endogenous retroviruses PERV-A/C from Göttingen minipigs.

BACKGROUND: Porcine endogenous retroviruses (PERVs) can infect human cells and pose a risk for xenotransplantation when pig cells, tissues or organs are transplanted to human recipients. Xenotransplantation holds great promise to overcome the shortage of human donor organs after solving the problems of rejection, functionality and virus safety. We recently described the transmission of a human-tropic recombinant PERV-A/C, designated PERV-F, from peripheral blood mononuclear cells (PBMCs) of a Göttingen Minipig (GöMP) to human 293 cells (Krüger et al., in Viruses 12(1):38, 2019). The goal of this study was to characterize PERV-F in more detail and to analyze the probability of virus isolation from other animals. METHODS: The recombination site in the envelope (env) gene, the long terminal repeats (LTR), the proteins and the morphology of the recombinant PERV-F were characterized by polymerase chain reaction (PCR), sequencing, Western blot analysis, immunofluorescence, and transmissible electron microscopy. Mitogen-stimulated PBMCs from 47 additional pigs, including 17 new GöMP, were co-cultured with highly susceptible human 293 T cells, and the PERV-A/C prevalence and PERV transmission was analyzed by PCR. RESULTS: PERV-F, isolated from a GöMP, is an infectious human-tropic PERV-A/C virus with a novel type of recombination in the env gene. The length of the LTR of PERV-F increased after passaging on human cells. In a few minipigs, but not in German landrace pigs, PERV-A/C were found. There was no transmission of human-tropic PERV-A/C from additional 47 pigs, including 17 GöMP, to human cells. CONCLUSION: These data show that human-tropic recombinant PERV-A/C proviruses can only be found in a very small number of minipigs, but not in other pigs, and that their isolation as infectious virus able to replicate on human cells is an extremely rare event, even when using highly susceptible 293 cells.

Human Endogenous Retroviruses as Gene Expression Regulators: Insights from Animal Models into Human Diseases.

The human genome contains many retroviral elements called human endogenous retroviruses (HERVs), resulting from the integration of retroviruses throughout evolution. HERVs once were considered inactive junk because they are not replication-competent, primarily localized in the heterochromatin, and silenced by methylation. But HERVs are now clearly shown to actively regulate gene expression in various physiological and pathological conditions such as developmental processes, immune regulation, cancers, autoimmune diseases, and neurological disorders. Recent studies report that HERVs are activated in patients suffering from coronavirus disease 2019 (COVID-19), the current pandemic caused by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) infection. In this review, we describe internal and external factors that influence HERV activities. We also present evidence showing the gene regulatory activity of HERV LTRs (long terminal repeats) in model organisms such as mice, rats, zebrafish, and invertebrate models of worms and flies. Finally, we discuss several molecular and cellular pathways involving various transcription factors and receptors, through which HERVs affect downstream cellular and physiological events such as epigenetic modifications, calcium influx, protein phosphorylation, and cytokine release. Understanding how HERVs participate in various physiological and pathological processes will help develop a strategy to generate effective therapeutic approaches targeting HERVs.

SARS-CoV-2 infection mediates differential expression of human endogenous retroviruses and long interspersed nuclear elements.

SARS-CoV-2 promotes an imbalanced host response that underlies the development and severity of COVID-19. Infections with viruses are known to modulate transposable elements (TEs), which can exert downstream effects by modulating host gene expression, innate immune sensing, or activities encoded by their protein products. We investigated the impact of SARS-CoV-2 infection on TE expression using RNA-Seq data from cell lines and from primary patient samples. Using a bioinformatics tool, Telescope, we showed that SARS-CoV-2 infection led to upregulation or downregulation of TE transcripts, a subset of which differed from cells infected with SARS, Middle East respiratory syndrome coronavirus (MERS-CoV or MERS), influenza A virus (IAV), respiratory syncytial virus (RSV), and human parainfluenza virus type 3 (HPIV3). Differential expression of key retroelements specifically identified distinct virus families, such as Coronaviridae, with unique retroelement expression subdividing viral species. Analysis of ChIP-Seq data showed that TEs differentially expressed in SARS-CoV-2 infection were enriched for binding sites for transcription factors involved in immune responses and for pioneer transcription factors. In samples from patients with COVID-19, there was significant TE overexpression in bronchoalveolar lavage fluid and downregulation in PBMCs. Thus, although the host gene transcriptome is altered by infection with SARS-CoV-2, the retrotranscriptome may contain the most distinctive features of the cellular response to SARS-CoV-2 infection.

Switching Sides: How Endogenous Retroviruses Protect Us from Viral Infections.

Long disregarded as junk DNA or genomic dark matter, endogenous retroviruses (ERVs) have turned out to represent important components of the antiviral immune response. These remnants of once-infectious retroviruses not only regulate cellular immune activation, but may even directly target invading viral pathogens. In this Gem, we summarize mechanisms by which retroviral fossils protect us from viral infections. One focus will be on recent advances in the role of ERVs as regulators of antiviral gene expression.

The Missing Expression Level-Evolutionary Rate Anticorrelation in Viruses Does Not Support Protein Function as a Main Constraint on Sequence Evolution.

One of the central goals in molecular evolutionary biology is to determine the sources of variation in the rate of sequence evolution among proteins. Gene expression level is widely accepted as the primary determinant of protein evolutionary rate, because it scales with the extent of selective constraints imposed on a protein, leading to the well-known negative correlation between expression level and protein evolutionary rate (the E-R anticorrelation). Selective constraints have been hypothesized to entail the maintenance of protein function, the avoidance of cytotoxicity caused by protein misfolding or nonspecific protein-protein interactions, or both. However, empirical tests evaluating the relative importance of these hypotheses remain scarce, likely due to the nontrivial difficulties in distinguishing the effect of a deleterious mutation on a protein's function versus its cytotoxicity. We realized that examining the sequence evolution of viral proteins could overcome this hurdle. It is because purifying selection against mutations in a viral protein that result in cytotoxicity per se is likely relaxed, whereas purifying selection against mutations that impair viral protein function persists. Multiple analyses of SARS-CoV-2 and nine other virus species revealed a complete absence of any E-R anticorrelation. As a control, the E-R anticorrelation does exist in human endogenous retroviruses where purifying selection against cytotoxicity is present. Taken together, these observations do not support the maintenance of protein function as the main constraint on protein sequence evolution in cellular organisms.

COVID-19 in Children: Expressions of Type I/II/III Interferons, TRIM28, SETDB1, and Endogenous Retroviruses in Mild and Severe Cases.

Children with the new coronavirus disease 2019 (COVID-19) have milder symptoms and a better prognosis than adult patients. Several investigations assessed type I, II, and III interferon (IFN) signatures in SARS-CoV-2 infected adults, however no data are available for pediatric patients. TRIM28 and SETDB1 regulate the transcription of multiple genes involved in the immune response as well as of human endogenous retroviruses (HERVs). Exogenous viral infections can trigger the activation of HERVs, which in turn can induce inflammatory and immune reactions. Despite the potential cross-talks between SARS-CoV-2 infection and TRIM28, SETDB1, and HERVs, information on their expressions in COVID-19 patients is lacking. We assessed, through a PCR real time Taqman amplification assay, the transcription levels of six IFN-I stimulated genes, IFN-II and three of its sensitive genes, three IFN-lIIs, as well as of TRIM28, SETDB1, pol genes of HERV-H, -K, and -W families, and of env genes of Syncytin (SYN)1, SYN2, and multiple sclerosis-associated retrovirus (MRSV) in peripheral blood from COVID-19 children and in control uninfected subjects. Higher expression levels of IFN-I and IFN-II inducible genes were observed in 36 COVID-19 children with mild or moderate disease as compared to uninfected controls, whereas their concentrations decreased in 17 children with severe disease and in 11 with multisystem inflammatory syndrome (MIS-C). Similar findings were found for the expression of TRIM-28, SETDB1, and every HERV gene. Positive correlations emerged between the transcriptional levels of type I and II IFNs, TRIM28, SETDB1, and HERVs in COVID-19 patients. IFN-III expressions were comparable in each group of subjects. This preserved induction of IFN-λs could contribute to the better control of the infection in children as compared to adults, in whom IFN-III deficiency has been reported. The upregulation of IFN-I, IFN-II, TRIM28, SETDB1, and HERVs in children with mild symptoms, their declines in severe cases or with MIS-C, and the positive correlations of their transcription in SARS-CoV-2-infected children suggest that they may play important roles in conditioning the evolution of the infection.

HIV-1 infection activates endogenous retroviral promoters regulating antiviral gene expression.

Although endogenous retroviruses (ERVs) are known to harbor cis-regulatory elements, their role in modulating cellular immune responses remains poorly understood. Using an RNA-seq approach, we show that several members of the ERV9 lineage, particularly LTR12C elements, are activated upon HIV-1 infection of primary CD4+ T cells. Intriguingly, HIV-1-induced ERVs harboring transcription start sites are primarily found in the vicinity of immunity genes. For example, HIV-1 infection activates LTR12C elements upstream of the interferon-inducible genes GBP2 and GBP5 that encode for broad-spectrum antiviral factors. Reporter assays demonstrated that these LTR12C elements drive gene expression in primary CD4+ T cells. In line with this, HIV-1 infection triggered the expression of a unique GBP2 transcript variant by activating a cryptic transcription start site within LTR12C. Furthermore, stimulation with HIV-1-induced cytokines increased GBP2 and GBP5 expression in human cells, but not in macaque cells that naturally lack the GBP5 gene and the LTR12C element upstream of GBP2. Finally, our findings suggest that GBP2 and GBP5 have already been active against ancient viral pathogens as they suppress the maturation of the extinct retrovirus HERV-K (HML-2). In summary, our findings uncover how human cells can exploit remnants of once-infectious retroviruses to regulate antiviral gene expression.